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Effects of resonant acoustic phonon scattering on magnetoresistivity are examined in two- 
dimensional electron systems at low temperatures by using a balance-equation magnetotransport 
scheme direct controlled by the current. The experimentally observed resonances in linear resistiv- 
ity are shown to result from the conventional bulk phonon modes in a GaAs-based system, without 
invoking leaky interface phonons. Due to quick heating of electrons, phonon resonances can be 
dramatically enhanced by a finite bias current. When the electron drift velocity increases to the 
speed of sound, additional and prominent phonon resonance peaks begin to emerge. As a result, re- 
markable resistance oscillation and negative differential resistivity can appear in nonlinear transport 
in a modest mobility sample at low temperatures, which is in agreement with recent experiments. 
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I. INTRODUCTION 

Low-temperature magnetoresistance oscillations re- 
lated to linear and nonlinear transport of electrons in 
high Landau levels of high- mobility two-dimensional sys- 



1,2,3,4,5,6,7,8,9 



or 



terns, induced by microwave radiation 
by direct current excitation ) 10 ' 11 ! 12 ' 13 ! 14 ' 15 ! 16 ! 17 have at- 
tracted a great deal of attention in the past few years. 
Despite the fact that detailed microscopic mechanisms 
are still under debate, there is almost no objection to 
refer these oscillatory magnetoresistances mainly to im- 
purity or disorder scatterings and the direct phonon con- 
tributions to resistivity are believed to be negligible in 
these systems at such low temperatures. 

Recently, the magnetophonon resonance in semicon- 
ductors, previously known to result from electron cou- 
pling with optic phonons and can be observed only at 
high temperatures and high magnetic fields , 18 ' 19 has 
been demonstrated to occur at temperatures as low 
as T ~ 3K and lower magnetic fields in GaAs-based 
heterosystemsj 20 ' 21 ' 22 ' 23 The resonant magnetoresistance 
was detected and referred to as electron scatterings by 
two leaky interface phonon modes that have sound ve- 
locities v s ~ 2.9km/s and v s ~ 4.4km/s ) 20 ' 21 ' 22 or by a 
single leaky interface phonon mode that has a velocity 
v s = 5.9km/si— Very recently, Zhang et al 2A found that 
acoustic phonon-induced resistance resonances are dra- 
matically enhanced in the nonlinear dc response and a 
finite current can strongly modify the phonon resonance 
behavior, transforming resistance maxima into minima 
and backi^i These phonon-related resistance oscillations 
remain poorly understood, especially the exact resonant 
condition, relative contributions of different modes and 
how they are affected by temperature and current. 

In this paper, we present a systematic analysis on non- 
linear magnetotransport in GaAs-based semiconductors 
with a microscopic balance-equation scheme directly con- 
trolled by the current, which takes into account electron 
couplings with impurity, bulk longitudinal and transverse 
acoustic phonons as well as polar optic phonons. Due to 



the quick rise of electron temperature, the phonon res- 
onances are dramatically enhanced by a finite current. 
When electron drift velocity v gets into the supersonic 
regime (v > v s , the speed of sound), additional mag- 
netophonon resonance peaks emerge. As a result, a re- 
markable resistance oscillation and a negative differential 
resistivity appear in the nonlinear magnetotransport in 
a modest mobility sample at low temperatures. 



II. BALANCE EQUATIONS FOR NONLINEAR 
MAGNETOTRANSPORT 



We consider a quasi-two-dimensional (2D) system con- 
sisting of N s electrons in a unit area of an x-y plane. 
These electrons, subjected to a uniform magnetic field 
B = (0, 0, B) along the z direction and a uniform elec- 
tric field E in the x-y plane, are scattered by random 
impurities and by phonons in the lattice. In terms of 
the center-of-mass momentum and coordinate defined as 

P = T,jPj\\ and R = N ^ J2j r j^ with Pj\\ = (Pjx,Pjy) 

and Tj — {xj , yj ) being the momentum and coordinate 
of the jth electron in the 2D plane, and the relative elec- 
tron momentum and coordinate p 1 ^ = Pj\\ — P/N s and 
r'j = Vj — R, the Hamiltonian H of this coupled electron- 
phonon system can be written as the sum of a center-of- 
mass part -ffcm, a relative electron part J/er ) 25 ' 26l27 ' 2S 



H cm = ——(P-N s eA{R)) — N s eE ■ R, (1) 
2iV s m 
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and a phonon part H p h — J2 q \ ^q\b f qX b qX , together 
with electron-impurity and electron-phonon interactions 
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as follows: 

#ci = Y. U ^h Z a)^ lqrra ^ qrR P qV (3) 
9||,o 

H cp = J2M(q,X)I(q z )(b q x + bl q x)^ R p qv (4) 

Here, -4.(t"), the in-plane component of the vector poten- 
tial of the uniform magnetic field, is linear in the spatial 
coordinate r = (x,y); m and m z are, respectively, the 
electron effective mass parallel and perpendicular to the 
plane; V{z) and V^(r^ — r'-, Zi, Zj) stand for the confined 

and Coulomb potentials, respectively p q ^ = ^■e l9|l " r J'H 
is the density operator of the two-dimensionalo (2D) rel- 
ative electrons; u(q^,z a ) is the effective potential of the 
ath impurity located at (r a , z a ) in the 2D Fourier space; 
b qX (b qX ) is the creation (annihilation) operator of the 
bulk phonon with wave vector q — (tjrii, q z ) — (q x , q y , q z ) 
in branch A that has an energy fi q \] M(q, A) is the ma- 
trix element of the electron-phonon interaction in the 
three-dimensional (3D) plane-wave representation; and 
I(q z ) is a form factor of the quasi-2D electron^! Here, for 
simplicity, we have assumed that the quasi-2D electrons 
occupy only the lowest subband and thus the subband 
summation indices in Eqs. ([3]) and ^ are neglected^ 

The separation of the electron Hamiltonian into a 
center-of-mass part and a relative electron part amounts 
to looking at electrons in a reference frame moving with 
their center of mass. The most important feature of this 
separation is that a spatially uniform electric field shows 
up only in 77 cm , and that H eT is the Hamiltonian of a 
many particle system subject to a perpendicular mag- 
netic field without the electric field. This enables us to 
deal with relative electrons in the magnetic field without 
tilting the Landau levels. The coupling between the cen- 
ter of mass and relative electrons is shown by the factor 
e ig H -R inside the momentum summation in H c \ and H cp . 
The moving center-of-mass assisted transitions of relative 
electrons between different Landau levels provide the ma- 
jor mechanism for the current-driven magnetotransport. 

Our treatment starts with the Heisenberg operator 
equations for the rate of change in the center-of-mass 
velocity V = —i[R, 77], and that for the relative electron 
energy 77 cr as follows: 

V = -i[V,H], (5) 
H et = -i[H ei ,H]. (6) 

When the electron-impurity and electron-phonon cou- 
plings are weak in comparison with the internal ther- 
malization of relative electrons and that of phonons, it is 
good enough to carry out the statistical average of the 
above operator equations to leading orders in H cl and 
77 cp . For this purpose, we only need to know the distri- 
bution of relative electrons and phonons without being 
perturbed by 77 oi or 77 op . The distribution function of 
the relative electron system described by Hamiltonian @ 
without an electric field should be an isotropic Fermi-type 



function with a single temperature T e . The phonon sys- 
tem, which is assumed to be in an equilibrium state, has 
a Bose distribution with lattice temperature T. Such a 
statistical average of the above operator equations yields 
the following force and energy balance equations in the 
steady state, which has a constant average drift velocity 
v: 

N s eE + N s e(v x B) + f(v) = 0, (7) 

v ■ f(v) +w(v) =0. (8) 

Here, f(v) = fi(v) + f p (v) is the frictional force experi- 
enced by the electron center of mass due to impurity and 
phonon scatterings, given by 

f i (v) = Y,Q\\\U(q\\)\ 2 n 2 (q h u; ), (9) 
911 

f p (v) = 2^<7|| \M(q,X)\ 2 \I(q z )\ 2 n 2 (q h Q qX +u ) 




and w{v) is the electron energy-loss rate to the lattice 
due to electron-phonon interactions with an expression 
obtained from the right-hand side of Eq. (|T0|) by replac- 
ing the qrii factor with fl q \. In these equations, ujq = q\\ -v, 
C/(<7||)| 2 is the effective average impurity scattering po- 
tential, and \M (q, X)\ 2 \I(q z )\ 2 is the effective coupling 
matrix element between a A-branch 3D phonon and a 
quasi-2D electron, /^(Qih is the imaginary part of the 
2D electron density correlation function at electron tem- 
perature T e in the presence of the magnetic field, and 
n(x) = l/(e x — 1) is the Bose function. The effect of 
interparticle Coulomb interaction is included in the den- 
sity correlation function to the degree of electron level 
broadening and screening. With the screening statically 
considered in the effective impurity and phonon poten- 
tials, the remaining LT2(q\\,Sl) function in Eqs. ([9]) and 
(fTU)) is that of a nonintcracting 2D electron gas in the 
magnetic field, which can be written in the Landau rep- 
resentation as2£ 

/Mini, fl) = ^^a,„'(^/2)/7 2 (n,n',/2), (11) 

" n,n' 

n 2 (n,n',fi) = ~ Jde [f(e) - f(s + Q)] 

x ImG„(e + Q) ImG n , (e), (12) 

where Ib = \J\j\eB\ is the magnetic length, C n ,n+i(Y) = 
nl[(n + /)!]- 1 r'e~ y [^(r)] 2 with L l n (Y) the associate 
Laguerre polynomial, /(e) = {exp[(e — p)/T c ] + 1} _1 
is the Fermi function at electron temperature T c , and 
ImG„(e) is the density-of-states of the broadened Landau 
level n. 

We model the electron density-of-states function with 
a Gaussian-type form for both overlapped and separated 



3 



Landau levels (e„ = nuj c is the center of the nth Landau 
level; ui c = eB/m is the cyclotron frequency) as follows: 29 

ImG„(e) = -(2^)5r- 1 exp[-2(e-e„) 2 /r 2 ]. (13) 

The half-width _T, or, the life-time or the quantum scat- 
tering time, t s — 1/2T, of the Landau level, which should 
be determined by electron-impurity, electron-phonon and 
electron-electron scatterings in the system, is magnetic- 
field B and temperature T dependent. We treat it as 
a semiempirical parameter, which will serve as the only 
adjustable parameter in the present investigation. 

At lattice temperature T, the energy balance [Eq. ([8])] 
yields the electron temperature T for a given carrier drift 
velocity v at a given magnetic field. Then, with this 
T e , the force balance [Eq. ([7])] determines the relation be- 
tween w, B and E, i.e., the longitudinal and transverse 
resistivities in the magnetotransport. 

Such a formulation indicates that the carrier drift ve- 
locity v is the basic physical quantity that controls the 
nonlinear magnetotransport. The frictional force f(v) 
and energy dissipation rate w(v) solely depend on the 
drift velocity v at a given magnetic field B, while the elec- 
tric field only plays a role in balancing the frictional force. 
Thus, the resistivity is directly determined by the scatter- 
ing mechanisms and by the drift velocity or the current 
density, rather than by the electric field. Equations ([7]) 
and ([5]) are conveniently applied to current-driven mag- 
netotransport of any configuration, in which the current 
is an experimentally directly controlled quantity. For an 
isotropic system where the frictional force is in the oppo- 
site direction of the drift velocity v and the magnitudes 
of both the frictional force and the energy-loss rate de- 
pend only on v = \v\, we can write f(v) = f(v)v/v 
and w(v) — w(v). In the Hall configuration with ve- 
locity v in the x direction [v — (v,0, 0)] or the current 
density J x = J = N s ev and J y = 0, Eq. ([7]) gives, af- 
ter T c [thus, the f(v) function] determined from Eq. ([5]), 
vf(v)+w(v) = 0, the transverse resistivity R yx = B/N s e, 
and the longitudinal resistivity R xx and the longitudinal 
differential resistivity r xx as 

R xx = -f(v)/(N 2 e 2 v), (14) 
r xx = -(df(v)/dv)/(Ny). (15) 

Note that, in principle, since the Landau-level broaden- 
ing and electron temperature are determined by the si- 
multaneous existence of all of the scattering mechanisms, 
contributions to the total resistivity from different scat- 
tering mechanisms are not independent. Nevertheless, 
it is still useful to formally write the total resistivity as 
a direct sum of separate scattering contributions by the 
respective components of the frictional force f(v). 

III. MAGNETOPHONON RESONANCE IN 
LINEAR MAGNETORESISTANCE 

In the numerical analysis we first concentrate on 
a GaAs-based heterosystem with carrier sheet density 



N s = 4.8 x 10 15 m~ 2 and zero-temperature linear mo- 
bility /io = 440 m 2 /Vs in the absence of a magnetic field, 
considering electron scatterings from bulk longitudinal 
acoustic (LA) phonons (one branch, via the deformation 
potential and piezoelectric couplings with electrons) and 
transverse acoustic (TA) phonons (two branches, via the 
piezoelectric coupling with electrons), as well as from re- 
mote and background impurities. The coupling matrix 
elements are taken to be well known expressions 27 with 
typical material parameters in bulk GaAs: electron effec- 
tive mass m — 0.067 m c (m is the free electron mass), 
longitudinal sound velocity v s i = 5.29 x 10 3 m/s, trans- 
verse sound velocity v st = 2.48 x 10 3 m/s, acoustic de- 
formation potential S — 8.5 eV, piezoelectric constant 
ei4 = 1.41 x 10 9 V/m, dielectric constant n = 12.9, 
and material mass density d = 5.31 g/cm 3 . We take a 
magnetic-field-dependent (B 1 ^ 2 ) Landau-level half-width 
as follows: 

r = [8aeu} c /Trmno(T)} 1/2 (16) 

which is expressed in terms of /j,q(T), the total linear 
mobility at lattice temperature T in the absence of the 
magnetic field, and a broadening parameter a to take 
into account the difference in the transport scattering 
time from the broadening-related quantum lifetime!^ 
A broadening parameter a — 1.5 is used in the calcu- 
lation in Sects. Ill and IV, which corresponds to taking 
Landau-level half- widths r — 1.8 and 2.4 K, respectively, 
at lattice temperatures T = 5 and 10 K for magnetic field 
B = 0.5T. 

The calculated total linear (v — > 0) magnetoresistiv- 
ity (differential resistivity) R xx = r xx is shown versus 
1/B at lattice temperatures T — 5K and T = 10 K 
in Figs. 1(a) and 1(b), together with separated contri- 
butions n from impurities, r p i from LA phonons, and r pt 
from TA phonons: r xx = r; + r p \ + r pt . The resistivity 
resonances clearly appear in r p i and r p t at both tem- 
peratures. This feature of magnetoresistivity stems from 
the property of 2D electron density correlation function 
772(911 1 ^q\)- In the case of low temperature (T e much less 
than the Fermi energy ep) and large Landau- level filling 
factor (v — ep/uJc ^> 1), 772(<7||, is a periodic function 
with respect to its frequency variable: n~2(q\\ , Q + lu> c ) — 
_/72 ( <? 1 1 , for any integer / of \l\ ^v. On the other hand, 
under the same conditions (low temperature T c <^ £p and 
large filling factor v 3> 1), the major contributions to the 
summation in Eq. (jlip come from terms n ~ n' ~ and 
then the function C n n > (x) has a sharp principal maxi- 
mum near x = 4t/i 30 i 31 Therefore, as a function of the in- 
plane momentum q\\, 7l2(g||, Q) sharply peaks at around 
g|| ~ 2fcp. In view of the existence of the form factor 
l-^fe)! 2 ! which is related to the wave function of the 
quasi-2D electron, only very small q z fcp) can con- 
tribute to the integration in Eq. (fTTJ)) . and thus, the con- 
tribution of the wave vector integration to the phonon- 
induced frictional force heavily weighs around q ~ 2fcp. 
As a result of this phase-space weight distribution and 
the 3D phonon dispersion [Q q \ = v s \q(X = 1, t)], when 



4 



: CO J CO 



1 



W = 4.8x1 ! 15 m" 2 m = 440 n Vs ( a ) 
v , = 5.29 km/s i/ = 2.48 km/s 








7= 10 K 


v=0 


(b) : 






r = r+r ,+r 

xx i pi F 

> ; 


4 




n 


\ X 








. . i . . . 


pt 


. . 1 . . 


2 



1 2 3 4 5 6 7 
S 1 (T 1 ) 

FIG. 1: (Color online) Linear (v — 0) resistivities (r mx , n, r p i, 
and r pt ) at (a) T = 5 K (a) and (b) T = 10 K for a GaAs-based 
heterosystem with carrier sheet density N B = 4.8 x 10 15 m -2 
and zero-temperature linear mobility fio = 440m 2 /Vs. The 
longitudinal and transverse sound velocities are taken to be 
v B i = 5.29 x 10 3 m/s and v st = 2.48 x 10 3 m/s. 
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FIG. 2: (Color online) (a) Differential resistivities r xx and 
(b) electron temperature T e are shown as functions of the 
magnetic-field strength B at lattice temperature T = 10 K 
under different values of the bias drift velocity v/vp. Here, 
vp is the Fermi velocity of the 2D electron system, and r xx 
curves of different 2v / v-p values are vertically offset for clarity. 



lu s \ = 2k-pv s \ = lu> c (I = 1,2,3,...), a quasi-2D relative 
electron in any Landau level can be resonantly scattered 
by absorbing or emitting a phonon and jumps across I 
Landau levels. The linear resistivity maxima then show 
up at 

e B A =w s \/w c = I = 1,2,3, ... (17) 

This is exactly what is seen in Figs. 1(a) and 1(b), as well 
as in the lower part of Fig. 2(a), where the r p \ maxima 
located at e s i = 1,2,3,4 and r p t maxima at e s % = 1,2 
are labeled. The total resistivity r xx peaks at around 
e s i = 1, 2, 3, 4, which is essentially determined by the res- 
onant scattering of LA phonons. These results are in 
good agreement with experiments ) 23 ' 24 which indicates 
that the observed low-temperature magnetophonon reso- 
nances in linear magnetoresistance are well explained by 
the ordinary single bulk LA phonon mode in GaAs with 
no need to invoke leaky interface modes. 

IV. MAGNETOPHONON RESONANCE IN 
NONLINEAR TRANSPORT 

A hnite current density J or a finite drift velocity v in 
the x direction has two major effects. First, it results in 
the electron heating and thus raises the rate of phonon 
emission. As a result, the phonon contributed resis- 
tivity and the oscillatory amplitude of magneto-phonon 



resonance is enhanced with increasing bias current den- 
sity. The electron temperature is determined by the fric- 
tional force f(v) and the energy-dissipation rate w(v) 
through the energy balance [Ea. ©] . Generally, longi- 
tudinal acoustic phonons give the dominant contribution 
to w(v) in GaAs-based systems at lattice temperatures 
considered in this paper (T < 10 K), and the polar op- 
tic (LO) phonons are usually frozen. However, electron 
scattering from LO phonons should still be taken into 
account when the bias current density becomes strong 
that electron temperature T c rises up to the order of 
20 K, at which a weak emission of LO phonons can take 
place. These emitted LO phonons, although giving lit- 
tle contribution to the resistivity itself, provide an ad- 
ditional efficient energy dissipation channel to prevent 
the continuing rise in electron temperature^ 2 - There- 
fore, we take the LO-phonon scattering (via the Frohlich 
coupling electrons, with an optical dielectric constant 
= 10.5) into account in the numerical calculation 
in nonlinear transport. The calculated electron tem- 
peratures T e in the case of T = 10 K at different bias 
drift velocities 2v/v F = 0.001,0.002,0.003,0.004,0.005, 
and 0.006 (vy is the Fermi velocity of the 2D elec- 
tron system), which correspond to current densities J = 
0.115, 0.231, 0.346, 0.462, 0.577 and 0.693 A/m, are shown 
in Fig. 2(b) versus the magnetic field B for the GaAs het- 
erosystem introduced above. 
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We see that, in addition to the concordant rise in elec- 
tron temperature with increasing current density, at each 
bias drift velocity, T c also exhibits resonance when chang- 
ing magnetic field. The oscillatory peak-valley struc- 
ture of T e periodically shows up, which has a period 
A(lo s \/lu c ) rs 1, indicating that it is also a magneto- 
phonon resonance mainly due to LA phonons. This T e 
resonance stems from the periodicity of f(v) and w{v) 
functions in the energy-balance equation (JS]) . Note that 
when the bias current density grows the oscillatory am- 
plitude of T c generally increases, but the positions of T c 
peaks (valleys) remain essentially the same in the B axis. 
By comparing with the magnetophonon resonant linear 
resistivity r xx [shown on the lower part of Fig. 2(a)], the 
oscillatory structure of T c peaks (valleys) appears roughly 
7r/ 2 phase shift down along the B axis. 

The differential resistivity r xx of the system at a given 
bias drift velocity v is obtained from the df(v)/dv func- 
tion through Eq. (|15[) . with the electron temperature T e 
determined above. Since the peak (valley) positions of T c 
are essentially fixed in the B axis, the variation of the re- 
sistivity maxima (minima) with changing bias velocity v 
is mainly determined by its direct effect in the df(v)/dv 
function, as reflected in the frequency shift ujq = qu-v in 
the argument of the electron density correlation function. 
Physically, because an extra energy qu • v is provided by 
the moving center of mass to the relative electrons during 
the scattering process, the transition rate of an electron 
from Landau level n to n' (n' can be equal to or not 
equal to n) experiences a change due to impurity and 
phonon scatterings, which shows up through the —qu ■ v 
shift in, for instance, the II2 (q\\ , O q \ — q\\ ■ v) function 
in — f p (i.e., the resistivity). The effect of such an en- 
ergy shift in the impurity scattering case has been shown 
to induce an oscillatory differential magnetoresistance n 
that is controlled by the following parameter : 10 i 12 i 14 

Sj = u>j/uJc, <jJj = 2k?v = y/Sw/N^J/e, (18) 

which has peak positions at around sj rs 0,1,2,.... In 
the case of acoustic phonon scatterings, the resistivity 
maxima are expected to occur near the possible integer 
values of parameter e s \ — £, , i.e., 



/ = 0,±1,±2, . 



(19) 



and this integer can be used as an identification for each 
magnetophonon resonance peak of the differential resis- 
tivity in nonlinear transport. When drift velocity v is 
smaller than the sound speed v s \ (ej < e s \), the res- 
onance condition (fT9"|) can be satisfied only for positive 
integers I = 1, 2, 3, 4, whence the effect of a small cur- 
rent J is roughly to change the resistivity maxima from 
e s \ Ri I in the linear case to 



(20) 
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This means that for a given index number I, the peak po- 
sition of magnetophonon resonance moves towards lower 



FIG. 3: (Color online) Differential resistivities (n, r p i, r pt , 
and r xx ) and electron temperature T e are shown as functions 
of Sj at fixed magnetic fields B = 0.27, 0.35, 0.47, and 0.72 T, 
respectively, at T — 5K. The number near the peak of r p i or 
r p t is the value of the integer I in the resonance condition (|19p 
of the corresponding differential resistivity peak. 



B with increasing current density, and for a given J the 
shifts of peak positions in the £ s a axis are larger for larger 
index I (lower B field) . These features are clearly seen in 
Fig. 2(a), where we show the calculated differential resis- 
tivity r xx as a function of the magnetic field at T = 10 K 
under different bias drift velocities from 2v/vp = to 
0.006 in steps of 0.0005, corresponding to current den- 
sities J = - 0.693 A/m in steps of 0.058 A/m. We 
see that though within certain v and -B-field ranges, e.g. 
2v/vp = 0.003 ~ .005 and B is around 6T, the peak 
positions of the resistivity may be somewhat influenced 
by the rapid change in the electron temperature because 
of the enhanced phonon contributions at higher T c ; the 
main trend of the resistivity peak shift with increas- 
ing J remains. The progress and movement of respec- 
tive peaks compare favorably with recent experimental 
observation.— 

When v becomes equal to or greater than the sound 
speed v s \ , the condition (fT9")) can be satisfied by I = and 
negative integers, which indicates the occurrence of ad- 
ditional magnetophonon resonance peaks under a strong 
dc excitation. Figure 3 presents the calculated differen- 
tial resistivity r xx and electron temperature T c at the 
lattice temperature T = 5 K as functions of the cur- 
rent density in terms of Sj = ujj/uj c at magnetic fields 
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FIG. 4: (Color online) Differential resistivity is plotted 
versus magnetic field B for fixed values of Sj = 0, 1, and 1.5 
at T = 5K. 



5 = 0.27,0.35,0.47, and 0.72 T. Due to the rapid rise 
in electron temperature T c in this fio = 440m 2 /Vs sys- 
tem, acoustic phonon contribution to resistivity becomes 
significant when Ej > 0.1. The prominent oscillations 
show up in resistivities r iy r p i, r pt , and r xx , all exhibiting 
a main oscillation period AEj ~ 1. However, since the 
peak positions of LA and TA phonon resistivities r p i and 
r pt in the Ej axis change with changing magnetic field, 
as indicated by condition (fT9| (sj ~ e s \ — I), the oscil- 
lating behavior of total resistivity r xx is strongly £?-field 
dependent. For instance, in the case of B = 0.35 T, the 
LA phonon resonance peaks 1,0,-1, and -2 are essentially 
in phase with the TA phonon resonance peaks 0,-1,-2 
and -3, and close to the maxima of the impurity resis- 
tivity n near Ej = 1,2,3,4, which lead to an enhanced 
r xx with maxima around these positions [Fig. 3(c)]. In 
the case of B = 0.47 T, where the LA resonance peaks 
1,0,-1, and -2 are around Ej = 0.5,1.5,2.5, and 3.5 
and the TA resonance peaks 0,-1,-2, and -3 are around 
Ej = 0.7, 1.7,2.7, and 3.7, the resulting r xx oscillation is 
out-of-phase with that of r; oscillation and has secondary 
extremes [Fig. 3(b)]. As a result, at fixed Ej, e.g., Ej = 2, 
r xx exhibits maxima for B = 0.35 and 0.72 T while min- 
ima for B = 0.27 and 0.47 T. Note that the resonant 
magnetophonon resistivity at v — v s \ (the I — peak) is 
generally the largest among all of the resistivity maxima 
in each type of phonon scattering, as can be seen from all 
four cases in Fig. 3. It comes from phonon-induced intra- 
Landau-level scatterings of electrons, which are allowed 
for all of the wavevectors, with the energy provided by 
the center of mass that has velocity v = v s \. 

At fixed Ej, phonon-induced resistivity is periodic in 
inverse B field with period Ae s \ » 1. The total resis- 
tivity r xx roughly follows this rule, as shown in Fig. 4, 
where r xx is plotted versus B for fixed Ej = 1 and 1.5. 
Both curves show remarkable oscillations that have the 
same periodicity but in- and out-of-phase as that of the 
zero-bias (sj = 0) case and with much enhanced am- 
plitudes. These features are just those observed in the 
experiments^ 
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FIG. 5: (Color online) Absolute R xx and differential r xx re- 
sistivities and their respective contributions from impurity, 
LA phonons, and TA phonons (Ri and n, R p \ and r p i and 
Rpt and r p t) and the electron temperature T e are shown ver- 
sus Ej at a lattice temperature T — 2K and a magnetic field 
B — 0.51 T for a GaAs-based heterosystem. The inset plots 
r xx at T = 2, 4, 6, and 8K with corresponding half-width of 
the Landau- level F. 



V. LOWER MOBILITY SYSTEMS 

Magnetophonon resonance can show up in resistivity 
at even lower lattice temperature as long as a finite cur- 
rent density is applied. In a system with lower mobility, 
the phonon contributions to resistivity can be greatly en- 
hanced by a modest finite current due to the rapid rise 
in the electron temperature. 

As another example, Fig. 4 plots the total absolute re- 
sistivity R xx and total differential resistivity r xx and their 
respective contributions from impurity, LA phonons, and 
TA phonons (Ri and R p \ and r p \ and R p t and r pt ), and 
the electron temperature T e , as functions of bias current 
density in terms of Ej at lattice temperature T — 2 K and 
magnetic field B = 0.51 T for a GaAs-based heterosys- 
tem that has carrier sheet density N s = 8 x 10 15 m -2 and 
zero-temperature linear mobility /io = 80m 2 /Vs in the 
absence of magnetic field. The Landau-level half-width 
r is taken to be 1.3 K (2T w 0.25w c at this magnetic 
field). 

Note that the phonon-induced absolute resistivities R p \ 
and R pt can become negative in this case. Although the 
total resistivity R xx remains always positive, the grow- 
ing negative r p \ and r pt with increasing current density at 
the initial stage greatly accelerate the drop in the total 
differential resistivity such that r xx goes down and be- 
comes negative at ej ~ 0.07, which is much smaller than 
it would be without phonon contribution (ej w 0.14), 
solely determined by the width of the Landau levelJ^ 
As a result, the total differential resistivity is negative 
within a wide current range 0.07 < ej < 0.42 before it 
becomes positive with further increasing current density. 
This may provide a possible explanation for the zero- 
differential resistance state recently observed at quite 
small dc bias in impure samples^ 
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As an important parameter in the present analysis, 
the half-width F of the Landau level, or the single par- 
ticle life-time or t s = 1/2 J 1 of the electron in the mag- 
netic field, which depends on all kinds of elastic and in- 
elastic scatterings induced by electron-impurity, electron- 
phonon and electron-electron interactions in the system, 
is temperature T dependent. At fixed B field, r certainly 
increases with increasing temperature due to enhanced 
electron-phonon and electron-electron scatterings. With 
this in consideration, the sharp drop in differential resis- 
tivity r xx at T — 2 K rapidly disappears when increasing 
temperature, as shown in the inset of Fig. 5, which is 
calculated with enlarged F at T = 4, 6, and 8 K. 



Good agreement between theoretical prediction and ex- 
perimental observation is obtained by using the well- 
defined bulk LA and TA phonon modes with no ad- 
justable parameter. We find that a finite dc excitation 
can not only greatly enhance the phonon resistivity by 
heating the electron gas, but can also induce additional 
and prominent phonon resonance peaks, which give rise 
to remarkable resistance oscillations and negative differ- 
ential resistivity in nonlinear transport. 



VI. SUMMARY 

In summary, we have examined magneto-phonon reso- 
nances in GaAs-based 2D systems at low temperatures. 
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tional Science Foundation of China and the Shanghai 
Municipal Commission of Science and Technology. 
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